Stimulated decelerating of diatomic molecules on multiple rovibrational transitions 

with coherent pulse trains 
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We propose a method of stimulated laser decelerating of diatomic molecules by counter- 
propagating 7r-trains of ultrashort laser pulses. The decelerating cycles occur on the rovibrational 
transitions inside the same ground electronic manifold, thus avoiding the common problem of radia- 
tive branching in Doppler cooling of molecules. By matching the frequency comb spectrum of the 
pulse trains to spectrum of the R-branch rovibrational transitions we show that stimulated cooling 
can be carried out on several rovibrational transitions simultaneously, thereby increasing number of 
cooled molecules. The exerted optical force does not rely on the decay rates in a system and can 
be orders of magnitude larger than the typical values of scattering force obtained in conventional 
Doppler laser cooling schemes. 
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While laser cooling is one of the key techniques 
of modern atomic physics (H-Q, neutral molecules are 
notoriously challenging to cool to ultracold tempera- 
tures. To accomplish this feat one must exercise pre- 
cise control over a multitude of internal degrees of free- 
dom. A breakthrough in molecular cooling techniques — 
even for diatomic molecules — is anticipated to enable 
substantial progress in quantum information process- 
ing Q , matter- wave interferometric sensors Q , quantum- 
controlled chemical reactions [6] and precision measure- 
ments @Q. 

To date, the coldest samples of diatomic molecules 
were obtained by assembling them from ultracold atoms 
via photo- or magneto- association. This approach pro- 
duced a gas of ground-state polar molecules near quan- 
tum degeneracy [10]. However, so far only bialkali 
molecules have been produced in this way, and the num- 
ber of molecules produced is fairly small (~ 10 4 ). Direct 
laser cooling of molecules could yield substantially larger 
samples for a wider range of species. 

In traditional Doppler cooling the radiative force orig- 
inates from momentum transfer to atoms from a laser 
field and subsequent spontaneous emission in random di- 
rections. Repeating this optical cycle tens of thousands 
of times can slow down thermal beams and cool atomic 
samples down to the Doppler limit (typically mK). Un- 
fortunately, most atoms and all molecules can radiatively 
decay to a multitude of states. Exciting population from 
all these lower-energy states requires a large number of 
lasers, which makes the conventional scheme impracti- 
cal. Only for a narrow class of molecules with highly- 
diagonal Franck-Condon overlaps (i.e., "near-cycling" 
electronic transitions), this branching problem can be 
mitigated [IlKlIj. 

Here we explore an alternative: employ absorption and 
stimulated emission on weak transitions inside the same 
molecular electronic manifold. The concept of "stim- 
ulated cooling" was put forward by Kazantsev [13| in 
1974; he proposed illuminating a two- level atom with a 



sequence of alternatively directed and oppositely detuned 
7r pulses. For an atom initially in the ground state, a 
7r-pulse impinging from the left would drive the popula- 
tion upward, while the 7r-pulse from the right cycles the 
population back to the ground state. The net change of 
atomic momentum is twice the photon recoil momentum 
p r = hk Cl where k c is the laser wave- vector. The advan- 
tage of this scheme over Doppler cooling is that the mo- 
mentum transfer can occur much faster than the radiative 
decay; this is crucial for the use of rovibrational transi- 
tions, for which Doppler cooling is impractical due to the 
very long radiative lifetimes. Stimulated cooling with 
counter-p ropaga ting 7r-pulses was explored in a number 
of works |14l4lq|. Closely related to the stimulated cool- 
ing by 7r-pulses is bichromatic force cooling (l7l-[2^ . Here 
the trains of Kazantsev's 7r-pulses effectively arise from 
the beating of two CW lasers of different frequencies. 
While the experiments so far have been limited to atoms, 
bichromatic force cooling of diatomic molecules on near- 
cycling electronic transitions was recently proposed [Hj. 

In this work we explore molecular slowing with coher- 
ent pulse trains. Pulse trains offer several advantages. 
First, the broadband structure of mode- locked lasers al- 
lows one to address multiple rotational levels simultane- 
ously. Secondly, it simplifies stabilization of the multiple 
frequencies involved, even for species for which building 
reference cells is impractical [26]. Finally, by manip- 
ulating the phases of successive pulse trains, as shown 
below, the capture range of velocities addressed can be 
tuned. Looking at the state-of-the-art frequency combs 
(FC) [I3-[3Q|, we note that a fiber-laser-based FC with 
10 W average power was demonstrated [27j and the au- 
thors argue that the technology is scalable above 10 kW 
average power. The spectral coverage was expanded from 
optical frequencies to ultraviolet and to mid-IR spectral 
regions o| . The high resolution quantum control via 
the combination of p ulse shaping and frequency comb 
was shown in [3ll434j. The experiments on line-by-line 
addressing have been done (32| . 



We focus on the transitions between the rovibrational 
levels (v, J) inside the ground electronic potential. The 
frequencies v V j of the i?-branch (J — >> J' = J + 1) 
electric-dipole-allowed transitions between rovibrational 
manifolds read 

~ + ( 3B «' " W + " B v) J 2 , (1) 

where is the frequency of the X(v, J = 0) — » 

X(v' , J' = 1) transition and B v are the rotational con- 
stants for the u-th vibrational level [Hj. This pattern 
matches the frequency comb (FC) spectrum which con- 
sists of a series of sharp equidistant peaks (teeth) located 
at 

Vn = Vc + ^rep X 71 - , (2) 

where n is an integer number, v c is the carrier frequency, 
the teeth-spacing z/ rep = 1/T is defined in terms of pulse 
repetition period T, and <E> is the carrier-envelope-offset 
(CEO) phase. We will focus on the v = Q —> v' = 1 tran- 
sitions. The molecular and the comb spectra could be 
matched by choosing v c = Jq \ and v rep = (3Bi — Bo)/n', 
with n' being an integer. The R-branch spectrum (pp) 
becomes Doppler-shifted for moving molecules; we will 
adjust $ to follow the Doppler shift of the band- head 
frequency. 

As an example, consider the LiCl molecule; here v c ~ 
19.29 THz with T = 0.95 ns for the number of teeth 
between nearby rovibrational transitions n' = 40. For 
the PbO molecule, v c « 21.63 THz with T = 0.98 ns for 
n' = 15. 

Based on these observations we propose the follow- 
ing scheme of stimulated cooling of molecules with pulse 
trains. As illustrated in Fig. HJa), the essential idea is 
to replace each of the Kazantsev's 7r-pulses by a train of 
N pulses, each of pulse area = ir/N. Interference of 
molecular probability amplitudes induced by the pulses 
leads to a frequency-dependence of the cooling force re- 
sembling the FC spectrum. Cooling will occur simul- 
taneously on several rovibrational transitions. The el- 
ementary cooling sequence (or cycle) will consist of two 
subsequent counter-propagating tt- trains. Each cycle will 
transfer twice the recoil momentum p r = hk c . To max- 
imize the optical force we assume that the time delay 
between trains is negligible. 

For a moving molecule, Doppler shifts due to the co- 
and counter-propagating trains have opposite signs. The 
counter-propagating/upward-stimulating train must be 
tuned below the resonance and the co-propagating train 
above the resonance. This can be attained by keeping the 
carrier frequencies of the two trains the same but having 
two different CEO phases. 

Before proceeding to the full analysis, note that the 
spectrum given in Eq. (pQ) is for the R-branch. To avoid 
driving P-branch transitions, one would need to elimi- 
nate the low-frequency (n < 0) half of the comb spectrum 
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FIG. 1: (Color online) (a) Stimulated cooling by 7r-trains. (b) 
R-branch rotational transitions form a spectrum that matches 
the frequency comb spectrum of pulse trains, (c) Reset- 
ing molecular population to the ground vibrational manifold 
X(v = 0, J) through the exited electronic state A with broad- 
band laser pulses. 

in Eq. (j2j), thus reducing the problem to a collection of 
(separate) two-level systems. This truncation may be ac- 
complished by installing a dispersive element and a static 
mask at the output of the cavity; the mask would cut out 
a predefined spectral window. Such a technique has been 
used by the Orsay group [36] in pulsed laser cooling of 
Cs2 molecules. 

Furthermore, we would like to avoid driving the v = 
1 — > v = 2 transitions, see Fig. QJb). Fortunately, 
due to anharmonicity, v\^\ < z^'j 1 , so the v\\ band- 
head is shifted to lower frequencies, already masked- 
out. Thereby only the higher- J rotational transitions will 
overlap with the truncated comb spectrum. The maxi- 
mum J of the v = 1 — >• 2 band overlapping with the 
comb spectrum could be estimated as J ma x = (^o'i — 
z^q'i)/(£>2 + B\) w u e x e /B e , where we used the conven- 
tional spectroscopic notation [35]. For LiCl and PbO, 
this limits us to J < 6 and J < 10 correspondingly. 
This limitation to low J is consistent with our neglect of 
higher-order non-rigid-rotor terms in Eq. (pQ). This limi- 
tation is of little consequence: at temperatures below 300 
K, the majority of the molecular population is in these 
low- J levels. Another peculiarity is that dipole matrix 
elements vary across rovibrational transitions, while the 
trains must execute the full tt rotation for all transitions 
simultaneously; to meet this requirement pulse shapers 
would need to be used so that the Rabi frequency of var- 
ious transitions remains the same. 

Formalism. - Since the FC fields couple only pairs of 
molecular levels, the problem is reduced to finding the 
time evolution and optical force on a two-level system. 
Due to the rapid nanosecond time-scale of 7r-train cycles 
we may neglect slow radiative decay of the vibrational 
levels (relevant lifetimes are on the order of a millisec- 
ond). We also approximate short non-overlapping laser 
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pulses by ^-functions. In this limit, the pulses are fully 
characterized by the pulse area 6 and phase <£. 

In the interaction picture, the propagator /7t ra in for 
a train of N pulses may be decomposed into a prod- 
uct of propagators due to individual pulses U p : t/train = 
Un "'U2U1. The pre-train (t = t$) and post-train 
(t = t J + NT) values of the density matrix are related 

as p(t£+NT) = Strain p(^o) ^train- Tne propagator 
across the p th pulse may be found analytically [37j : U p = 
I cos (6/2) + icr p sin (6/2), with cr p = cr x cosr]L,R(t p ) — 
cr y sin r}L,R(t p ), where a XjV are the Pauli matrices. Here 
t p is the arrival time of the p-th pulse and L (left) and R 
(right) label the co- and counter-propagating 7r-trains in 
Fig. [I][a). The phase rj LiR (t) = StT k c z(t) - <$> L ' R (t) is 
the cumulative phase of the laser field experienced by the 
moving molecule, z(t) « vt. Focusing on a target velocity 
vo, we may redefine § L ' R (t) = (5 =p k c v (t))t + ^ R (t), 
where ^ ,R (t) is the control phase that we will tune to 
optimize the cooling process. Then 



VL,n(t) = Tk c (v - v (t))t - ^ R (t) . 



(3) 



At this point we have a prescription for evolving the 
density matrix over time. The last needed ingredient is 
the expression for the mechanical momentum transferred 
to the two-level system. The fractional momentum kick 
due to a single train may be expressed in terms of the 
excited state population difference at the end and at the 
beginning of the train -Ap train /p r = p ee (t£ + NT) - 

Pee {to) [23. 

For the elementary cycle 
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Pee (t+ + 2NT) + Pee (to) - 2 Pee (t+ 



NT) . 

(4) 
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FIG. 2: Dependence of the fractional momentum kick per 
cooling cycle on the Doppler defect phase. Here N c = 500 
and 7r-trains are composed of four 7r/4 pulses. The overview 
plot in panel (a) was generated while the velocities were kept 
constant and control phases ,R = 0. Panel (b) zooms in 
onto the central tooth of panel (a). 

Based on the described formalism we developed a 
Mathematica package to simulate stimulated cooling by 
coherent pulse trains. We start by discussing our com- 
puted dependence of momentum kick per cycle on the 
Doppler "defect" phase k c (v—vo(t))T, shown in Fig.[2fa). 



First of all, Ap is a periodic function, reflecting the un- 
derlying periodicity of the FC spectrum. The transferred 
momentum spikes at values of the phase which are mul- 
tiples of 2tt. At these points the maximum transferred 
momentum is limited to twice the recoil momentum, as 
expected from Eq. (J2J). This happens only when during 
the cycle the system starts and ends in the ground state, 
with a full transfer of the population to the excited state 
by the first train, i.e., when the 7r-train conditions are 
satisfied. Qualitatively, at such values of the phase the 
probability amplitudes transferred by subsequent pulses 
interfere constructively. 

The detailed profile of the momentum transfer shown 
in Fig. [2fb) depends on the number of pulses inside tt- 
trains and the number of cycles. For a fixed number 
of pulses the calculated profile displays a very compli- 
cated substructure resulting from intricate interferences 
of probability amplitudes driven by the multitude of 
pulses. In accord with the time- frequency uncertainty 
principle, as the number of pulses grows larger, the peaks 
become narrower effectively reducing the capture veloc- 
ity range to zero. Because of this effect, it is natural to 
wonder how to increase the number of cooled molecules 
at the end of the process. This goal can be attained by 
imposing phase relation between the last and the first 
pulses of two subsequent tt- trains: 



$f' L (W,l) = -$c'^M,iv) +^ 



L,R( 



(5) 



Here tM,n is the arrival time of the n-th pulse of the M-th 
train. The origin of the phase relation (|5]) may be readily 
understood using the Bloch sphere visualization method 
as illustrated in (supplemental Material at [URL will be 
inserted by publisher]. As the molecules slow down, the 
phases have to track (t) as proposed in Ref. [37| ; a typ- 
ical rate of phase-tuning is d^/dt ~ p^/(hm)/N, where 
m is the molecular mass. Experimentally the required 
phase-tuning can be attained with electro-optical modu- 
lators. 

In Fig. [3ja) we investigate the time evolution of the 
velocity distribution for LiCl molecules. We suppose that 
initial cryogenic molecular beam has a 140 m/s forward 
velocity and a rotational temperature of « IK [38]. The 
forward velocity distribution is modeled as a Gaussian 
with half- width equal to 20 m/s. 

Initially all molecules are in the ground rovibrational 
manifold with (J < J max) and they enter the interac- 
tion volume at the same moment of time (t = 0). The 
laser phase follows the time evolution of vq, so the maxi- 
mum of the optical force is always at vo(t). To maximize 
the velocity capture range, the control phase changes as 
Eq. (j5j). The strongly peaked force (Fig. [2j) leads to the 
compression of velocity distribution near its maximum. 
However since the force profile narrows with increasing 
number of cycles, the width of the maximum of veloc- 
ity distribution starts to decrease too. The area under 
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FIG. 3: (Color online) Time evolution of velocity distribu- 
tion of the ensemble of decelerated diatomic molecules LiCL 
All molecules are initially thermally distributed over the ro- 
tational states of ground vibrational state of the ground elec- 
tronic state. Initial forward velocity distribution of cryogenic 
molecular beam is modeled by the gaussian with center at 
140 m/s and half- width equal to 20 m/s. Single pulse area 
— 7r/4, pulse repetition period T = 1 ns. (a) Velocity 
distributions for the stimulated cooling scheme with Doppler 
phase tracking and train-to-train phase correction, (b) The 
same as in (a) but with additional population reset applied 
every 10,000 cycles. 



the maximum also decreases, that means that the effec- 
tive number of cooled molecules grows smaller. This is 
caused by a build-up of destructive interferences at the 
wings of the force profile. 

To improve the final number of decelerated molecules 
one could simply stop the trains and wait until all 
the molecules radiatively decay into the ground ro- 
vibrational manifold and then restart the cooling pro- 
cess. This would clear out all the unwanted superpo- 
sitions. However the required time-scales are too long. 
Instead, one could employ the much faster radiative de- 
cay of the upper electronic states, Fig. Die). In this 
process additional broad-band laser pulses would drive 
the transitions: X(v,J) — >• A(v',J'), with the lower- 
frequency transitions X(v = 0, J < J m ax) — >• A(v' ', J') fil- 



tered out. Due to the spontaneous decay of the A(v' , J') 
levels, most of the population after several absorption- 
radiative decay cycles gravitates towards the manifold 
X(v — 0, J < Jmax) used as a starting point for stim- 
ulated cooling. We simulated such a reset /stimulated 
cooling scheme in Fig. [3f b). Here the stimulated cooling 
process is the same as in Fig. [3ja), with the addition of 
resetting pulses applied every 10,000 cycles. We observe 
that while the peak of velocity distribution becomes nar- 
rower over time, its height increases so that the area (i.e., 
the number of molecules) increases. 

So far we neglected fine (or hyperfine) structure of 
molecular states. These could complicate the analysis as 
the pulses could couple several levels. Still one could iso- 
late two-level transitions using chiral pulse trains. For ex- 
ample, for the singlet ground state diatomic molecule one 
could use the <r + (<j~ ) circularly-polarized trains. Then a 
subset of two- level transitions (v = 0, J, F = 1+ J, Mp = 
I+J)±+(v = 1, J' = J+l, F' = J+ J+l, M F , = I+J+l) 
becomes isolated. 
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